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Faust Signal Processing Libraries 

• oscillator. lib — signal sources 

• filter.lib — general-purpose digital filters 

• effect.lib — digital audio effects 
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Highlights of Additions Since LAC-08 
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• oscillator.lib 

o Filter-Based Sinusoid Generators 
o Alias-Suppressed Classic Waveform Generators 

• filter.lib 

o Ladder/Lattice Digital Filters 
o Audio Filter Banks 


• effect.lib 

o Biquad-Based Moog VCFs 
o Phasing/Flanging/Compression 
o Artificial Reverberation 
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Moog Voltage Controlled Filters (VCF) 



Overview 


effect.lib 

• Moog VCF 

• phasing/flanging 

• reverberation 

filter.lib 


oscillator.lib 
Conclusion 


• moog_vcf _2b = ideal Moog VCF transfer function factored 
into second-order “biquad” sections 

o Static frequency response is more accurate than 
moog_vcf (which has an unwanted one-sample delay in 
its feedback path) 

o Coefficient formulas are more complex when one or both 
parameters are varied 

• moog_vcf _2bn = same but using normalized ladder biquads 

o Super-robust to time-varying resonant-frequency 
changes (no pops!) 

o See Faust example vcf_wah_pedals.dsp 

















Moog VCF 


Moog VCF 

moog.vcf(res,fr) 


moog_vcf _2b (res, f r) 
moog vcf 2bn(res, f r) 


See Faust example vcf _wah_pedals. dsp 

analog-form Moog VCF 

res = corner-resonance amount [0-1] 

f r = corner-resonance frequency in Hz 

Moog VCF implemented as two biquads (tf 2) 

two protected, normalized-ladder biquads (tf 2np) 
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Phasing and Flanging 


Phasing and Flanging 

vibrato2jnono(...) 
phaser2_mono (. . .) 
phaser2_stereo(. . .) 

flangerjnono(...) 
flanger stereo(. . .) 


See Faust example phaser_flanger .dsp 
modulated allpass-chain (see effect. lib for usage) 
phasing based on 2nd-order allpasses (see effect. lib 1 
stereo phaser based on 2nd-order allpass chains 
mono flanger 
stereo flanger 

















Artificial Reverberation (effect. lib) 



Overview 


effect.lib 

• Moog VCF 

• phasing/flanging 

• reverberation 

filter.lib 


oscillator.lib 
Conclusion 


• General Feedback Delay Network (FDN) Reverberation 

See Faust example reverb .designer .dsp 

• Zita-Revl Reverb (FDN+Schroeder) by Fons Adriaensen 
(ported to Faust) 

See Faust example zita_revl. dsp 
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Overview 


effect.lib 

filter.lib 

• ladder/lattice 

• normalized ladder 

• filter banks 

oscillator.lib 
Conclusion 


Ladder/Lattice Digital Filters (filter.lib) 

• Ladder and lattice digital filters have superior numerical 
properties 

• Arbitrary Order (thanks to pattern matching in Faust) 

• Arbitrary (Stable) Poles and Zeros 

• All Four Major Types: 


o Kelly-Lochbaum Ladder Filter 
o One-Multiply Lattice Filter 
o Two-Multiply Lattice Filter 
o Normalized Ladder Filter 
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Normalized Ladder Digital Filters (filter. lib) 



Overview 


effect.lib 

filter.lib 

• ladder/lattice 

• normalized ladder 

• filter banks 


oscillator.lib 


Conclusion 


Advantages of the Normalized Ladder Filter Structure: 

• Signal Power Invariant wrt Coefficient Variation 
Extreme Modulation is Safe 

• Super-Solid Biquad (sweep it as fast as you want!): 


tf2snp() 


“transfer function, 2nd-order, s-plane, normalized, protected” 
• See Faust example vcf_wah_pedals.dsp 
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Ladder and Lattice Digital Filters 


Lattice/Ladder Filters 

iir_lat2(bcoeffs,acoeffs) 
iir_kl(bcoeffs,acoeffs) 
iir_latl(bcoeffs,acoeffs) 
iir_nl(bcoeffs,acoeffs) 
tf2np(b0,bl,b2,al,a2) 

nlf2(f,r) 


two-multiply lattice digital filter 
Kelly-Lochbaum ladder digital filter 
one-multiply lattice digital filter 
normalized ladder digital filter 
biquad based on stabilized 
second-order normalized ladder filter 
second-order normalized ladder digital filter 
special API 
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Overview 


effect.lib 

import( 

filter.lib 

• ladder/lattice 

bcoeffs 

• normalized ladder 

• filter banks 

acoeffs 

oscillator.lib 

Conclusion 

process 


filter.lib"); 

= (1,2,3); 

= ( 0 . 1 , 0 . 2 ); 

= impulse <: 

iir(bcoeffs,acoeffs), 
iir_lat2(bcoeffs,acoeffs), 
iir_kl(bcoeffs,acoeffs), 
iir_latl(bcoeffs,acoeffs) 


: > 
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Audio Filter Banks (filter. lib) 



Overview 


effect.lib 

filter.lib 

• ladder/lattice 

• normalized ladder 

• filter banks 

oscillator.lib 
Conclusion 


• “Analyzer” = Power-Complementary Band-Division 
(e.g., for Spectral Display) 

See Faust example spectral_level.dsp 

• “Filterbank = Allpass-Complementary Band-Division 
(Bands Summable Without Notch Formation) 

See Faust example graphic_eq.dsp 


Julius Smith 


• Filterbanks in filter. lib are implemented as analyzers in 
cascade with delay equalizers that convert the 
(power-complementary) analyzer to an 
(allpass-complementary) filter bank 
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oscillator.lib 
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oscillator.lib 


Reference implementations of elementary signal generators: 


effect.lib 


filter.lib 

oscillator.lib 

• sinusoids 

• oscb 

• oscr 

• oscs 

• oscw 

• virtual analog 

• sawN 

• sawtooth examples 

• pink noise 

Conclusion 


• sinusoids (filter-based) 

• sawtooth (bandlimited) 

o pulse-train = saw minus delayed saw 
o square = 50% duty-cycle pulse-train 
o triangle = (leakily) integrated square 
o impulse-train = differentiated saw 
o (all alias-suppressed) 


• pink-noise (1//noise) 
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Overview 


effect.lib 
filter.lib 

oscillator.lib 

• sinusoids 

• oscb 

• oscr 

• oscs 

• oscw 

• virtual analog 

• sawN 

• sawtooth examples 

• pink noise 

Conclusion 


Sinusoid Generators in oscillator. lib 


oscb 

“biquad” two-pole filter section 
(impulse response) 

oscr 

2D vector rotation 
(second-order normalized ladder) 
provides sine and cosine outputs 

oscrs 

sine output of oscr 

oscrc 

cosine output of oscr 

oscs 

state variable osc., cosine output 
(modified coupled form resonator) 

oscw 

digital waveguide oscillator 

oscws 

sine output of oscw 

oscwc 

cosine output of oscw 
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Block Diagrams 


Inspect the following test program: 


effect.lib 


filter.lib 


oscillator.lib 

• sinusoids 

• oscb 

• oscr 

• oscs 

• oscw 

• virtual analog 

• sawN 

• sawtooth examples 

• pink noise 

Conclusion 


import("oscillator.lib"); 

freq = 100; 

process = oscb(freq), 

oscrs(freq), 
oscs(freq), 
oscws(freq); 
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• sinusoids 

• oscb 
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• oscs 
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• virtual analog 
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Sinusoidal Oscillator oscb 

oscb (impulsed direct-form biquad) 

• One multiply and two adds per sample of output 

• Amplitude varies strongly with frequency 

• Numerically poor toward f req=0 (“dc”) 

• Nice choice for high, fixed frequencies 
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Sinusoidal Oscillator oscr 
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filter.lib 
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• sinusoids 

• oscb 

• oscr 

• oscs 

• oscw 

• virtual analog 

• sawN 

• sawtooth examples 

• pink noise 

Conclusion 


oscr (2D vector rotation) 

• Four multiplies and two adds per sample 

• Amplitude is invariant wrt frequency 

• Good down to do 

• In-phase (cosine) and phase-quadrature (sine) outputs 

• Amplitude drifts over long durations at most frequencies 
(coefficients are roundings of s = sin(2*PI*freq/SR) 
and c = cos(2*PI*freq/SR), so S 2 + c V 1) 

• Nice for rapidly varying frequencies 

















Sinusoidal Oscillator oscs 



Overview 


effect.lib 
filter.lib 

oscillator.lib 

• sinusoids 

• oscb 

• oscr 

• oscs 

• oscw 

• virtual analog 

• sawN 

• sawtooth examples 

• pink noise 

Conclusion 


oscs (digitized “state variable filter”) 

• “Magic Circle Algorithm” in computer graphics 

• Two multiplies and two additions per output sample 

• Amplitude varies much less with frequency than oscr 

• Good down to do 

• No long-term amplitude drift 

• In-phase and quadrature components available at low 
frequencies (exact at dc) 


• Nice lower-cost replacement for oscr when amplitude can 
vary slightly with frequency and exact phase-quadrature 
outputs are not needed 
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Sinusoidal Oscillator oscw 
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effect.lib 
filter.lib 

oscillator.lib 

• sinusoids 

• oscb 

• oscr 

• oscs 

• oscw 

• virtual analog 

• sawN 

• sawtooth examples 

• pink noise 

Conclusion 


oscw (2nd-order digital waveguide oscillator) 

• One multiply and three additions per sample (fixed frequency) 

• Two multiplies and three additions when frequency is changing 

• Same good properties as oscr, except 

o No long-term amplitude drift 

o Numerical difficulty below 10 Hz or so (not for LFOs) 

o One of the two state variables is not normalized (higher 
dynamic range) 


• Nice lower-cost replacement for oscr when state-variable 
dynamic range can be accommodated (e.g., in VLSI) 

















Virtual Analog Waveforms in oscillator. lib 
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effect.lib 


filter.lib 


oscillator.lib 

• sinusoids 

• oscb 

• oscr 

• oscs 

• oscw 

• virtual analog 

• sawN 

• sawtooth examples 

• pink noise 

Conclusion 


imptrain(freq) 
squarewave(freq) 
sawtooth(freq) 
sawN(N,freq) 


periodic impulse train 
zero-mean square wave 
alias-suppressed sawtooth 
order N anti-aliased saw 


• sawtooth and sawN based on “Differentiated Polynomial 
Waveform” (DPW) method for aliasing suppression 

• sawN uses a differentiated polynomial of order N 
Increase N to reduce aliasing further 

• Default case is sawtooth = saw2 = sawN(2) 

(sounds quite good already!) 

• Bandlimited square, triangle, and pulse-train derived as linear 
filterings of bandlimited sawtooth 
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Faust Source for sawN 


sawN(N,freq) = sawl : poly(N) : D(N-l) : gate(N-l) 
with { 

pOn = float(ml. SR) /float(freq); // period in samples 

lfsawpos = (_,l:fmod) ~ +(1.0/p0n); // sawtooth in [0,1) 

sawl = 2*lfsawpos - 1; // zero-mean, amplitude +/- 1 

poly(l,x) = x; poly(2,x) = x*x; 

poly(3,x) = x*x*x - x; ... 

diffl(x) = (x - x’)/(2.0/p0n); 

diff(N) = seq(n,N,diff1); // N diffls in series 

D(0) = 

D(1) = diffl/2.0; 

D(2) = diff(2)/6.0; 

• • • 

gate(N) = *(1@(N)); // blanks startup glitch 

>; 
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Sawtooth Examples 


Faust Examples Using Bandlimited Sawtooth saw2 

(saw2(freq) = sawl(freq) <: * <: -(mem) : *(0.25 J *SR/freq) ;) 

• <f aust>/examples/graphic_eq. dsp 

• <faust>/examples/gate_compressor.dsp 

• <faust>/examples/parametric_eq. dsp 

• <f aust>/examples/phaser_f langer. dsp 

• <f aust>/examples/vcf _wah_pedals. dsp 
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Pink Noise 


• Pink noise has the same power in every octave, making it perceptually 
more uniform than white noise 

• oscillator. lib implements pink_noise (“1 // noise”) 
(approximately) as white noise through a three-pole, three-zero MR filter 
that approximates a 1// power response: 


pink_noise = noise : 

iir((0.049922035, -0.095993537, 0.050612699, -0.004408786), 

(-2.494956002, 2.017265875, -0.522189400)); 

• This filter was designed using invf reqz in Octave (matlab) by fitting 
three poles and zeros to a minimum-phase I/a/ 7 amplitude response 
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Overview 


effect.lib 


filter.lib 


oscillator.lib 


• Main developments in Faust signal-processing libraries 
oscillator | filter I effect. lib since LAC-08 were 
summarized 


Conclusion 

• Conclusion 

• Acknowledgments 


• Ongoing goal is accumulation of reference implementations in 
music/audio signal processing 
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